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VII. QUANTUM MECHANICAL LANGEVIN EQUATIONS 38 
A Rudimentary Resevoir Problem ; 3940 

Consider a complete system which consists of a single simple harmonic oscillator (the 
system component of "interest") coupled to a reservoir of many simple harmonic 
oscillators (the system component treated statistically) 



"The System" 
A single simple 
harmonic oscillator 



^nt 



eraction 




Reservoir 
of 

many harmonic 
oscillators 



In particular, the complete interacting system has a simple Hamiltonian of the form 

M = M laser + ^ res+ ^ int 

= h Q da+Y^h (Oj b f jbj +^ti [gj a< b j + g* b] a} 



[VII-1] 



where a is the system variable of interest (the "state of the system") and the 's are the 

other system variables (the "reservoir states"). System operators commute with reservoir 
operators at given time, we may easily generate the following set of equations of motion: 



Much of what follows draws heavily on material in the Arizona Books — i.e. 

1. M. Sargent III, M. O. Scully and W. E. Lamb, Jr., Laser Physics, Addison- Wesley (1974) 

2. P. Meystre and M. Sargent III, Elements of Quantum Optics, Springer- Verlag (1992) 

3. Weng W. Chow, Stephan W. Koch and Murray Sargent III, Semiconductor-Laser Physics, 
Springer- Verlag (1994) 

This model is particular valuable in treating electromagnetic intermode coupling problems, but it also provides 
guidence in treating more general problems which incorporate atomic states. 
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d(t) = i [?{,a(t)]/h= -z'Q a(t)-iY,gj bj(t) [ VII-2a ] 

j 

bj(t) = i [tt,b j (t)]/h= -/to. bj(t) - ig* a(t) [ VII-2b ] 
Formally integrating Equation [ VII- 2b ], we obtain 41 

bj( t) = bj ( t Q ) expj - i (i> j (t - f „)] - i g* J dt a{t) expj - i (Oj(t- /)] [ VII-3 ] 

The first term in this equation describes the free evolution of the reservoir states in the 
absence of any interaction with the system and the second gives the modification of this 
free evolution as a consequence of the coupling to the system. Inserting Equation [ VII-3 ] 
into Equation [ VII-2a ] we obtain 

d(t) = -iQ.a(t)-Y} gj f J dt'a (j')exp[-i 1 CO . (t- t)] 

1 [ VII- 

- i L gj b M ex p[-' a j { f - r o)] 

Here the second summation describes a source of fluctuations arising from the free 
evolution of the reservoir states and the first is a feedback through the reservoir which 
might be described as a radiation reaction. To remove the rapid variation in the system 
variable we transform to the Heisenberg interaction picture - i. e. we take 42 

a(t) = A(t) exp[-z Q f] [ VII-5 ] 

As is pointed out in Section 14-3 of P. Meystre and M. Sargent III, Elements of Quantum Optics, Springer- 
Verlag (1992), the statement that "System operators commute with reservoir operators at given time. .." is a 
very critical point. From Equation [ II-3 ], we see that as time evolves, the reservoir operators acquire some of 
the characteristics of the system operator. Thus, it is crucial in this development to maintain a strict and 
consistent ordering of operators at different times! 

The transformation does not alter the system's communtation relationship, so, of necessity, [A{t),A\t)]=l. 
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sothat M*)=-Yi 8j f J ' d/A(/)exrj[-/(co.-n)(*-^] + F(f) [VII-6] 

where F(t) = -i£ 8j bj(t ) exp[i (Q-co .) ? )] [ VII-7 ] 

j 

is the so called noise operator. 

The first term on the RHS of Equation [ VII-6 ] may be reordered as 
Y}gj | 2 J ' dfAtfad-ifa-Ot) (t- 0] 

j 

= \' dt'A (t) £| gj | 2 exrj-/ (to. -a) (f- 0] 



[ VII-8a ] 



By assumption, the frequencies of the reservoir oscillators are closely spaced and widely 

distribution so that the function L| gj | exp[-z (co^ -Q) (/- f')] is very sharply peaked at 

j 

t'= t and has a width of the order of the inverse of the reservoir's frequency bandwidth. 
Thus, A{t') may taken out of the integration and the limits of the integration may be 

extended to infinite ~ i. e. 



Yjgj f J ' ^'A(f')exp[-/(co,-Q) (t- 0] 

j 

^ A (') LI gj f |^^exp[-i (co, -a) (t- t)] 



[ VII-8b ] 



Again following arguments explicated by Heitler, 43 we see that 



43 In Chapter II, Section 8 of W. Heitler, The Quantum Theory of Radiation (3rd edition) 
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J dt exp[-/(a-0) .) (t- tj\ =>jc 8(a-a) .)- Tr —!— [ VII-9 ] 

If we neglect the imaginary part of Equation [ VII-9 ], which leads to a Lamb-type 
frequency shift, Equation [ VII-6 ] leads directly to a quantum mechanical version of the 
classical Langevin equation 44 — viz. 

A(t) = -^A(t) + F(t) [VII-10] 

where y= 2tc p (ft) \ g(Q) f . [ VII- 1 1 ] 

is the quantum mechanical Langevin drift coefficient. In optical applications the average 
intensity or number operator — i.e. (a'(?)A(^ — is the experimental quantity of 

greatest interest. A useful equation of motion for that operator may obtained by making use 
the Langevin equation ~ i.e. Equation [ VII-10 ]. To that end we first write 

j- t [A\t) A (t))= -y(A\t) A(t)) + {F\t) A (t)) + (A\t) F(t)) [ VII-12 ] 
To obtain (F\t) A (t)) and (A\t) F (t)) we make use of the identity 



44 In the classical Langevin theory of Browian motion, the equation of motion for drift velocity of a particle 
suspended in a liquid is given by 

mil (t) =—m Yu (t) + F u (t) 
where Y is a damping constant and F v (t) is a random, rapidly fluctuation force with zero mean value. In this 
classical problem the famous fluctuation-dissipation theorem is expressed as 

Y=^D uu =^\ + ~dx{F u (0)F u (T)) 
where D uu is the so called Fokker-Planck diffusion coefficient. 

R. Victor Jones, May 2, 2000 76 



The Interaction of Radiation and Matter: Quantum Theory Page 77 

A(t)=A(t-At)+\' dt'AfA [VII-13] 

so that 

lF*(t) A (t)) = (F\t) A it- At))+ I' dt (F\t) A it)) 

•><-& i [VII-14] 

= [F\t)A{t- At))- I \'_ A dt'(F\t)A(t)) + \jt'[F{t)F{t')) 

The first term on the RHS vanishes since cause cannot precede effect. That is, A(t) 
cannot be correlated with a future value of the fluctuating force. Similarly, in the second 
term ^F ] (t) A(t)) vanishes except in the infinitesimal small inteval where t' = t. 

Therefore 

(F\t)A(t)) = J ' Jt'(F\t) F(t)) . [ VII-15a ] 

If the random force is stationary in time 

(F\t) A (t)) = J ° a dx (F\t) F(t+x)} [ VII-15b ] 

and, finally, if A? is long compared to the random force correlation time 

(F\t) A (t)) = J ° dx (F\t) F (t+x)) = ^ J + °°dx (F\0) F (x)) . [ VII-15c ] 

Therefore Equation [ VII- 12 ], the equation of motion for the average number operator, 
becomes 

| (A\t) A(t))=-y (A\t)A(t))+ J Jx (f\0),F(x)) [ VII-16 ] 
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which is one form of the famous fluctuation-dissipation theorem. 

If the reservoir is in thermal equilibrium, we may directly evaluate the noise operator 
correlation to wit. 

( F VM0L rm = I 8* & (bp) b k (0)\ herm exp[-/(n-co,) t] exp[/(Q-(D t ) t] [ VII- 

i,k 

Since the operators of the reservoir commute we can write 

(bmbM) =(/»,) 8 >t 

\ ]\ I k\ 'I therm » 1 'therm ]K 

9 f„ 1 
= ™ In^exjj-n. (3 ft ©J }• [ VII-18 ] 

dp U J 

= {exp[(3^Q) J ]-l}" 1 6 J , 

and Equation [ VII- 17 ] becomes 

(*M*(0L.=Ekf («,),», «iHn-»,)(«-0] 

= 2ti I g(a) | 2 p (a) (n(Q)) ;/ _5( ? -0 

= Y(^)L rm 5(^0 [VII-19] 

\ I therm V 7 

/ \ Y 

where ( D t ) = — ( n (Q) ) is of the form of a diffusion coefficient in the 

\ / therm 4 therm 

classical Langevin theory. Two-time correlation functions of this form are characteristic of 
Markoffian random processes and represent fluctuations in a reservoir with essentially 
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zero memory. Substitution of this Markoffian correlation function into Equation [ VII- 16 ] 
yields an Einstein relation of the form 



General Resevoir Problem 

As the previous section demonstrates, the effects of random fluctuations in an assemblage 
of harmonic oscillators may be accounted for in precise detail. The problem can be 
analyzed completely since the first term on the RHS of Equation [ VII-6 ] — the damping 
term — contains only the system variable. In general, this feedback term will included 
reservoir variables as well and the simplified model breaks down. The general case may, 
however, be analyzed if we take the simpler case as a guide. In particular, let us suppose 
that a set of systems operators { \{t) } are coupled to the reservoir and satisfy a set 

Langevin equations of motion as 



2l d t ) = — 1 

\ A * A I therm rff^ 



/ th, 



:erm 



[ VII-20 ] 



[ VII-21 ] 



where D^(t) is the drift term and F^t) is the Markoffian noise operator appropriate to the 
system operator A^(7) • From the earlier results, we are, thus, assuming that the two-time 

correlation of the random fluctuating force is of the form 



(WV(0)=2^ V 5(^-0 [VII-22] 



Again using the identity 



\(t) = \(t-At) + \ t Jt'A ll (t') 



[ VII-23 ] 



we find 
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{\(t)Fit))=(\(t-At)FXt)) 

c, r, [VII-24] 

+ \^dt'(DMF y (t)) + \^dt'(FMF y (t)) 

As discussed in connection with Equation [ VII- 14 ], the first and second term on the RHS 
in Equation [ VII-24 ] vanishes since cause cannot precede effect. Therefore 

(\{t)Fit))=\ljt{F^t)FXt)) [VII-25] 

and by Equation [ VII-22 ] 

[\{t)K{t))=D^ [VII-26] 
We may then form the following equation of motion 

jJ\\ A v) = (\ A >A\ A \) [vn27] 
= ( D , A ,) + ( F , A v ) + (\D v ) + {A ll F v ) 

which leads to a generalized Einstein relationship 

2D ^ = J t \\ A v)- (D, A v ) - (\ D v ) [ VII-28 ] 

With this equation it is possible to calculate the diffusion coefficients from the drift 
coefficients, provided one can independently calculate the equation of motion for 'A^ A v ^. 

We can also use Equation [ VII-21 ] to calculate the expectation values of two-time 
correlation functions. In paryicular, if multiply it by A v (? / ) where t < t and average we 

find 
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j t (\(t)Mt))=(DM 4M>+teM 4,(0) [ vn-29 ] 

Since in the Markoff approximation, the present system operator A v (/) cannot know about 
the future noise source F^(t) , (F^t) A^ty) vanishes and 

j t M*W)) = to to 4(0) [vn-30] 

which is the so called quantum regression theorm which states that the two-time correlation 
function ^ (?) A v (tf} satisfies the same equation of motion as the single-time (A^ [t)\ . 
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Semiconductor Langevin Equations 

As discussed earlier, the effective Hamiltonian for optical interactions in a semiconductor 
may be written as 



M eft- - 



Thus, the equation of motion motion for the effective dipole operator b _ k & k is given by 



\L n 2 e 


t 


h 2 k 2 ^ 













d ir -i 

6 r ar=-\!H eff ,b r a r 

dt ^ * 



=-(y+ ico-J 6 _ i a- k +ig- k a(al a- k + 6 ^6 + f 6 ^ 



[ VII-30 ] 



h 2 k 2 h 2 k 2 

where y is the damping, C0- = £ H — H — , and ^ „ is the k-dependent noise 

2m 2m, ~ k k 

ec hv 

operator. We transform the various operators to a frame rotating at the electromagnetic 
field frequency (0 r — viz. 



6 -i fl k = ° k exp[-co r ?] 
a(t)=A(t) exp[-co r t] 
fe.m= fot= F o* ex P[-£0, '] 



[ VII-31a] 
[VII-31b] 
[VII-31c] 



and find the transformed equation of motion in the rotating frame — i.e. the Heisenberg 
interaction picture ~ as 

j f o- =- ( y+ / CO- - / tO r ) o- + ig - A (a la- k + 6 t 6 _- - 1) + [ VII-32 ] 
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Applying the generalized Einstein relation ~ i.e. Equation [ VII-28 ] ~ to the operator 



we obtain 



Similarly, 



2fl» fl =(y+i(o- k -i(» r )o k 33 +(y- /co- + /(o r ) < + a 



« 2y a 



33 



1D a at = (Y+ i®*- i »r) <4 +(Y- * » k + °00 + 



[ VII-33a ] 



[ VII-33b ] 



~2ya* 

where aj, = (l - a { - 6 _ f - 6 _-) . 

The equation of motion for the laser field operator in a frame rotating at the laser frequency 
(£> r may be written as 



v \ r co , , 



4'HEft<*i(') + *W [VII-34] 



where Q. is the cold cavity resonant frequency and where we have replace the damping 
parameter y in equation [ VII- 10 ] by the experimental parameter (o/Q . If we consider the 
Equation [ VII-32 ] in the quasi- equilbrium limit, we find 



=[Y+ *e>i- i(O r ]~ l [ig- k A{t)(ala- k + B \6 _--l) + F 0k (?) 
= 2)(co- - to r ; y) \jg - A{ t) - a k QQ ) + F 0| ( *)] 



[ VII-35 ] 
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where Z((o- k -(o r ;y) is the complex Lorentzian demoninator. 45 Thus in the quasi- 
equilbrium limit, Equation [ VII-34 ] becomes 



r 



^(^)=-^ + ^-«J-L|d 2 ^ ro k- w ^Y)(o3 3 -o^)jA(?) + J F A (?) [VII-36] 



where 



[ VII-37 ] 



Following arguments similar to those used in deriving Equation [ VII- 16 ], we easily derive 
an equation of motion for the laser photon number operator ~ viz. 



jMX'MO)- 



V ' i J [VII-38a] 



+ \~dx(F}{0),F A {%)) 



or 



-{ A \t)A(t)) = -\^-(R 33 -R 0Q ) 



Too 



(A\t) A(t))+ \*~dx (Fl(0),F A (x)) I VII-38b ] 



where 



• k 



[ VII-39a ] 



45 It may be recalled that in the notes entitled Lasers: Models and Theories for convenience we indroduced the 
Lorentzian functions -- viz. 

( T{u — v; w) = [z [u — v) + w] m e complex Lorentzian denominator 

£[u - v; w) = w 2 [(w - v) +w 2 | the dimensionless Lorentzian function. 

R. Victor Jones, May 2, 2000 84 



The Interaction of Radiation and Matter: Quantum Theory Page 85 

and R QQ = ^Lkf^K - * ^) °oo ■ [ VII-39b ] 



Y 



k 



Assuming that the cavity damping and carrier scattering reservoirs are uncorrelated, the 
two-time correlation function may be written as 

(Ft(>)F A (>)) = (FX>)F(>)) 

Y) J/K-» r ;l) (j>M F„.(0) 1 VII - 4 ° a 1 

k k' 

and if the noise operators for different k are uncorrelated 

W m) = F{ t ))A lh^k-^ Y ) ( F : k (t) F Ok .(0) . [ vn-40b ] 

' k 

Finally, using Equations [ VII- 19 ] and [ VII-33a ] we obtain 

(#(*K(0)= IW^L^+^Eld'^-^'^J^jsC'-O • [vii-40c] 

and, thus, the "A A" diffusion coefficient for the laser field is 

2^ = -{«ML rm +*33 • [VII-41a] 
By a similar argument, the " A A " diffusion coefficient is 

2D a* tim +l ] + R n- [VII-41b] 



Using Equation [ VII-41a ], Equation [ VII-38b ] becomes 
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CO 



— - (7? 33 - R Q() ) 



33 



[ VII-42a ] 



Similarly 



r co , k 

~ [ R 33 ~ R 00) 



(A(t)A\t)) + ^[[n(co)) 



+l]+R m 

therm J uu 



[ VII-42b ] 



As check on the consistency of these result, we see that the form of Equations [ VII-42a ] 
and [ VII-42b ] guarantees that the commutator [ A(t), A l (t)] is independent of time - viz. 



dt 



GO 



— - (i? 33 - R QQ ) 



([A(t),A\t)])+^ + (R QQ -R 33 )= [ VII-43 ] 



to incorporate saturation effects into the noise spectrum we need to write a Langevin 
equation for the carrier-density operator - viz. 



— fl i! fl k = — ["^ eff ' fl k fl k]~'~ {imping, pumping, and noise terms} 



dt 



= A] [l- a( aj - T- a{ a k b \ 6 _ k - y nr a{ a- k 



[ VII-44 ] 



a l a - k -(4 fl J J - 1 k a l 6 - k a ~ £ fl k ^-J +F e 



where A k e [l - a £ a -] i s the pumping due injected carriers, T- a I a- b _ + - b - is the radiative 
recombination (spontaneous emission) term, Y„ r fl t fl k is the nonradiative recombination 



term, r « 



is the relaxation due carrier-carrier scattering, and is the 
carrier noise operator. In the quasi-equilibrium approximation — i.e. from Equation 
[ VII-35 ] — we may sum Equation 111-44 ] to obtain a Langevin equation for the total 
carrier-density operator 
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— N = A- — Yll a^-a-b^-b - - y nr N 
dt V - 

2A 1 A 2 - - [VII-45a] 

-— TT-Lkf 4 a) k- o) r ;Y)(°33-o!; )+F iV 

' k 

where 

F N (*)= Fht) + ^[gtA\t) ^a)--co r ; Y ) F 0i (t) - g- k A(t) <D*(<o- k -<o r ;y) F^t)] [VII-45b] 

V k 

Noise Spectra: 

To calculate noise spectra we make use of semiclassical transformations - viz. 

A{t) = (e^/eJ exp[-/<K?)] [ VII-46a ] 

A\t) = (E{t)/<£ a ) exp[i<|)(f)] [ VII-46b ] 

E[t) = 'E ar ^A\t)A{t) [VII-46C] 
^t) = (i/2) ln[A( t)/A\ t)] [ VII-46d ] 

where 2^ = [h (0 r /e V] is the electric field per photon. Using the expression 

d' = Yd [vn-47] 

po £ -34 3a; LV114/J 

we see that 
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D EE- D A A Ja +L>AA " ^7 ~ 



a A 3 A 1 



4E Z 



AA* D A , A +A'AD AA ]~^[D A , A + D AA ] 



w 3A 1 3A dA 3A 1 



1 1 £ 2 r i 

= TAA 7 Da ^ + 4A T A ~4~fu DAtA + jC>AAt -l 



[ VII-48a ] 



[ VII-48a ] 
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